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[t]  This  paper  presents  an  analysis  of  the  wave  number  spectra  of  intermediate-scale 
waves  (wavelengths  between  0.02  and  6  m)  covering  a  wide  range  of  wind  and 
background  wave  conditions  under  various  sea-state  conditions.  The  main  result  of  the 
analysis  is  that  the  dependence  of  the  dimensionless  wave  spectrum  on  the  dimensionless 
wind  friction  velocity  follows  a  power  law  function.  The  coefficient  and  exponent  of  the 
power  law  function  vary  systematically  with  the  wave  number.  The  wave  number 
dependence  of  the  coefficient  and  exponent  serves  as  an  empirical  parameterization  for 
computing  the  wave  number  spectra  of  intermediate-scale  waves  at  different  wind  speeds. 

Calculation  of  the  mean  square  slope  from  the  resulting  wave  number  spectrum 
confirms  that  intennediate-scale  waves  are  the  dominant  contributor  of  the  ocean  surface 
roughness.  A  simple  formula  is  presented  for  calculating  the  band-pass  filtered  mean 
square  slope  of  the  ocean  surface  for  remote  sensing  applications. 

Citation:  Hwang,  P.  A.  (2005),  Wave  number  spectrum  and  mean  square  slope  of  intermediate-scale  ocean  surface  waves, 
J.  Geophys.  Res..  110,  C10029,  doi.10.1029/2005JC003002. 


1.  Introduction 

[2]  Surface  waves  a  few  centimeters  to  a  few  meters  long 
(referred  to  as  the  intermediate-scale  waves  hereafter)  are 
generally  believed  to  be  die  dominant  contributor  of  the 
ocean  surface  roughness.  Owing  to  the  large  Doppler 
frequency  shift  when  traditional  point-measurement  systems 
are  employed,  the  measurement  accuracy  of  intermediate- 
scale  waves  in  the  ocean  is  difficult  to  assess.  So  far,  reliable 
measurements  of  intennediate-scale  waves  in  the  ocean 
environment  are  essentially  nonexistent.  As  a  result,  verifi¬ 
cation  of  the  assumption  regarding  intermediate-scale  waves 
and  the  ocean  surface  roughness  based  on  field  data  is 
not  yet  presented.  Most  ocean  wave  data  focus  on  the 
energy-containing  frequency  band  in  the  neighborhood  of 
the  spectral  peak  of  the  sea  surface  displacement.  A  small 
fraction  of  the  wave  data  focuses  on  the  length  scale  of 
gravity-capillary  waves,  using  mainly  optical  scanning 
sensors  to  conduct  spatial  measurements  to  avoid  the  tricky 
problems  associated  with  the  Doppler  frequency  shift  [e.g., 
Hwang  et  al.,  1996;  Hara  et  al.,  1998],  Extending  the 
optical  scanning  techniques  to  wave  components  several 
meters  long  for  field  applications  is  prohibitively  expensive 
and  impractical.  Hwang  et  al.  [1996]  suggested  that  the 
problem  of  Doppler  frequency  shift  can  be  alleviated  by 
using  a  free-drifting  measurement  technique.  The  effective¬ 
ness  of  removing  Doppler  frequency  shift  by  free-drifting 
operations  was  demonstrated  by  comparing  the  wave 
number- frequency  spectra  measured  by  linear  arrays  of  wire 
gauges  deployed  in  free-drifting  and  fixed-station  config¬ 
urations  [Hwang  and  Wang,  2004]  and  reliable  conversion 
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of  the  wave  spectrum  from  frequency  to  wave  number 
domain  can  be  extended  into  intermediate-scale  wave  com¬ 
ponents.  A  more  quantitative  discussion  of  the  Doppler 
frequency  shift  is  given  in  Appendix  A. 

[3]  Hwang  and  Wang  [2004]  presented  an  analysis  of 
field  measurements  of  intermediate-scale  waves.  Following 
the  suggestion  by  Phillips  [1984],  the  data  were  processed 
to  investigate  the  empirical  functional  dependence  of  the 
dimensionless  spectrum,  B(k),  on  the  dimensionless  wind 
friction  velocity,  u*/c,  where  k  is  tire  wave  number  and  c  is 
the  wave  phase  speed.  The  result  shows  that  B(k)  can  be 
expressed  as  a  power  law  function  of  m*/c.  The  coefficient 
and  exponent  of  the  power  law  function  are  wave  number 
dependent.  In  this  paper,  using  the  above  result  serving  as 
an  empirical  parameterization  function,  the  wave  number 
spectra  of  intennediate-scale  waves  at  different  wind  speeds 
are  quantified  (section  2).  The  properties  of  the  mean  square 
slopes  integrated  over  different  wave  number  ranges  are 
investigated  in  section  3.  A  summaiy  is  given  in  section  4. 

2.  Empirical  Parameterization 

[4]  In  a  discussion  of  the  source  function  balance  of  short 
ocean  surface  waves  that  are  important  to  remote  sensing 
applications,  Phillips  [1984]  emphasized  that  the  knowledge 
on  the  variation  of  the  wave  number  spectrum  with  wind 
speed  can  provide  valuable  information  on  the  properties  of 
the  dissipation  function.  Hwang  and  Wang  [2004]  reported 
an  analysis  of  the  wave  spectra  collected  in  the  ocean  using 
a  free-drifting  technique.  The  spectra  are  divided  into  two 
different  groups:  wind  seas  and  mixed  seas.  The  range  of 
wind  speeds  is  3.6  to  14.2  m/s  for  the  former  group  (291 
cases),  and  2.6  to  10.2  m/s  for  the  latter  group  (106  cases). 
For  each  wave  component,  the  scatterplot  of  B(k)  as  a 
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Figure  1.  (a)-(c)  Examples  of  the  scatterplot  of  B(u*/c)  for  a  given  wave  number  component,  showing 
a  power  law  dependence.  The  wave  number  is  marked  at  the  upper  left  comer  of  each  panel.  The  size  of 
the  plotting  symbol  is  proportional  to  the  data  density  at  the  plotting  coordinates,  (d)  A0  and  (e)  a0  of  the 
power  law  function  derived  from  least  squares  fitting  of  field  data.  See  color  version  of  this  figure  in  the 
HTML. 


function  of  u*/c  shows  a  general  power  law  dependence,  as 
illustrated  by  the  examples  shown  in  Figures  la-lc, 

(it*  \  /«* 

*(-;*)  =4>(*)(t)  •  (i) 

The  coefficient  and  exponent  of  the  power  law  function  for 
each  wave  component  can  be  calculated  from  least  squares 
fitting  of  the  ocean  data.  The  variations  of  A0  and  a0  with  k 
are  depicted  in  Figures  Id  and  le,  respectively.  Hwang  and 
Wang  [2004]  presented  fifth-order  polynomial  fitting 
functions  for  Aq  and  a0  but  the  polynomial  functions 
obviously  missed  many  fine  features  of  the  variations  (see 
their  Figure  3).  Here  the  computation  is  earned  out  by  using 
a  lookup  table  interpolating  through  the  measured  data 
(Table  1),  the  interpolated  curves  capture  most  of  the  subtle 
details  of  the  variations  in  A0  and  a0  (compare  the 
continuous  curves  and  discrete  data  in  Figures  ld-le). 
For  the  range  1  <  k  <  2  rad/m,  the  scatter  of  the  mixed-sea 
data  is  large  and  the  wind-sea  data  in  this  range  are  used  for 
both  wind  seas  and  mixed  seas.  The  wave  number  spacing 
shown  in  Table  1  is  logarithmic.  Other  interpolation 
methods  such  as  linear  and  cubic-spline  interpolations  in 
the  logarithmic  or  linear  spacing  were  tested  and  they 
produced  only  small  differences  in  the  final  construction  of 
A0  and  n0  curves  or  the  subsequent  calculation  of  B(k)  at 
different  wind  speeds. 

[5]  Figure  2  displays  the  computed  wave  number  spectra 
over  a  range  of  wind  speeds  using  the  parameterization 


Table  1.  Lookup  Table  of  Aa  and  a0  for  Wind  Seas  and  Mixed 
Seas 


k,  rad/m 

Wind  Seas 

Mixed  Seas 

103.4„ 

«0 

lOUfn 

CIO 

1.00 

37.32 

1.04 

37.32 

1.04 

1.26 

18.60 

0.75 

18.60 

0.75 

1.58 

16.28 

0.67 

16.28 

0.67 

2.00 

14.92 

0.62 

14.92 

0.62 

2.51 

13.09 

0.56 

13.09 

0.56 

3.16 

10.30 

0.45 

10.11 

0.46 

3.98 

8.97 

0.35 

9.98 

0.43 

5.01 

8.21 

0.24 

10.06 

0.36 

6.31 

10.27 

0.32 

10.52 

0.31 

7.94 

9.83 

0.32 

9.08 

0.24 

10.00 

10.71 

0.33 

9.70 

0.23 

12.59 

11.99 

0.37 

10.06 

0.22 

15.85 

12.31 

0.37 

11.50 

0.28 

19.95 

12.24 

0.36 

10.72 

0.24 

25.12 

12.49 

0.41 

10.60 

0.22 

31.62 

12.70 

0.47 

11.63 

0.33 

39.81 

12.38 

0.49 

11.59 

0.36 

50.12 

11.80 

0.52 

10.80 

0.37 

63.10 

11.29 

0.59 

10.42 

0.43 

79.43 

10.82 

0.69 

9.77 

0.49 

100.00 

9.77 

0.76 

9.24 

0.59 

125.89 

8.88 

0.87 

8.13 

0.63 

158.49 

7.58 

0.98 

7.08 

0.71 

199.53 

6.31 

1.10 

5.88 

0.78 

251.19 

4.85 

1.24 

4.36 

0.86 

316.23 

3.61 

1.39 

3.34 

0.62 
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Figure  2.  B(k)  computed  from  parameterization  function  (1)  using  table  lookup  for  A0  and  ci0  (Table  1): 
(a)  Wind  seas  and  (b)  mixed  seas  plotted  in  semilogarithmic  scales  (area  conservative  representation  of 
the  spectrum  of  mean  square  slope),  (c)  The  normalized  difference  of  the  saturation  spectra  in  wind  seas 
(denoted  by  subscript  w )  and  mixed  seas  (denoted  by  subscript  s).  (d)  B(k),  with  the  kink  near  k  =  1  rad/m 
smoothed,  presented  in  log-log  scales.  Wind  seas  are  plotted  with  curves  and  mixed  seas  with  markers. 
See  color  version  of  this  figure  in  the  HTML. 


function  (1)  for  the  wave  number  range  1  <k<  316  rad/m. 
The  semilogarithmic  plots  in  the  top  panels  represent  the 
area-conservation  representation  of  the  mean  square  slope 
of  the  ocean  surface;  that  is,  the  area  under  the  curve  over  a 
range  of  wave  numbers  represents  the  mean  square  slope  of 
the  corresponding  spectral  wave  components,  because  for 
the  omni-directional  spectrum  B(k)  =  kSt(k),  where  S\  is  the 
omni-directional  wave  slope  spectrum  (Appendix  of  Hwang 
et  al.  [1993]).  The  results  for  the  wind  seas  are  plotted  in 
Figure  2a  and  those  for  the  mixed  seas  in  Figure  2b.  The 
normalized  spectral  difference  between  the  mixed  seas  and 
wind  seas  is  shown  in  Figure  2c.  In  mixed  seas,  the 
background  ocean  waves  modify  the  spectrum  of  interme¬ 
diate-scale  waves  in  a  somewhat  complex  way.  For  waves 
shorter  than  about  1  nr,  tire  background  waves  enhance  the 
spectral  densities  of  intermediate-scale  waves  at  lower  wind 
speeds  but  reduce  their  spectral  densities  at  higher  wind 
speeds.  For  longer  intermediate-scale  waves,  the  effect  of 
background  waves  is  opposite  to  that  for  the  shorter  waves. 
The  magnitude  of  the  spectral  fluctuations  deviating  from 
the  wind-sea  condition  is  usually  within  about  20  percent 
for  most  intermediate-scale  wave  components  in  the  avail¬ 
able  data.  The  critical  wind  speed  separating  these  two 
opposite  trends  is  somewhat  higher  than  6  m/s  (Figure  2c). 
Interestingly,  the  wind  speed  7  m/s  is  frequently  associated 
with  the  inception  of  more-intensive  breaking  events  in  the 
ocean  and  the  transition  of  the  surface  roughness  condition 


from  smooth  or  transitional  to  hydrodynamically  rough.  The 
results  shown  in  Figure  2c  may  reflect  the  influences  of 
breaking  wave  and  surface  roughness  conditions  on  the 
intermediate-scale  waves.  The  exact  effects  of  wave  break¬ 
ing  and  roughness  conditions  on  the  properties  of  short-  and 
intermediate-scale  waves  in  the  ocean  are  not  well  under¬ 
stood.  Theoretical  investigations  on  related  subjects  have 
indicated  that  the  mechanisms  of  wave-turbulence  interac¬ 
tion  and  wave  breaking  exert  strong  impact  on  the  air-sea 
momentum  transfer  and  wind-wave  generation  [e.g.,  Makin 
and  Kudryavtsev,  1999,  2002],  From  the  empirical  data 
gathered  here,  it  is  found  that  the  response  of  wave 
components  between  1  and  2  m  long  to  the  background 
waves  is  opposite  to  that  of  shorter  waves.  The  data  quality 
of  wave  components  longer  than  3m(l</c<2  rad/m)  in 
the  present  dataset  is  not  very  good  for  the  mixed-sea 
conditions,  judging  from  the  large  data  scatter  of  A0  and 
a0  in  this  wave  number  range  shown  in  Figures  ld-le,  and 
it  is  difficult  to  determine  the  influence  of  background 
waves  in  mixed  seas  on  the  longer  components  of  interme¬ 
diate-scale  waves. 

[s]  For  the  wave  number  range  kp<k<  1  rad/m,  where  kp 
is  the  wave  number  at  the  peak  of  the  surface  displacement 
spectrum,  the  equilibrium  spectral  function  is  assumed, 

Se{k)==  bu*g~05k~2'5  =  b—lc\  (2) 
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Figure  3.  Mean  square  slopes  integrated  from  kp  to  different  upper-bound  wave  numbers,  ks.  Field  data 
from  sun  glitter  and  altimeter  scattering  analyses  are  also  presented  for  comparison.  See  color  version  of 
this  figure  in  the  HTML. 


where  the  spectr  al  constant  b  «  5.2  x  I  O'2  is  used,  and  g  is 
the  gravitational  acceleration.  The  corresponding  dimen¬ 
sionless  spectrum  is 

Be{k)=b(^>).  (3) 

[7]  The  results  shown  in  Figure  2  illustrate  that  inter¬ 
mediate-scale  waves  are  the  dominant  contributor  of  the 
ocean  surface  mean  square  slope.  Numerical  integration 
shows  that  more  than  77  percent  of  the  mean  square  slope 
is  contributed  by  waves  between  0.02  and  6  m  for  wind 
speeds  less  than  20  m/s.  The  peak  of  the  mean  square 
slope  spectrum  represented  by  B(k )  is  near  20  rad/m  at 
3  m/s,  and  moves  toward  higher  wave  number  as  wind 
speed  increases,  reaching  to  about  200  rad/m  at  18  m/s. 
The  dropoff  of  the  spectral  densities  at  the  higher  wave 
number  end  (waves  shorter  than  about  2  cm)  is  quite 
steep.  There  is  an  apparent  kink  at  k  =  1  rad/m  in  merging 
B(k)  and  Ba(k),  possibly  due  to  the  relatively  simple 
equilibrium  spectral  model  used  here.  This  is  not  consid¬ 
ered  to  be  a  significant  drawback  for  the  discussion  of  the 
ocean  surface  roughness  because  the  contribution  of  long 
gravity  waves  to  the  overall  mean  square  slope  of  the 
ocean  surface  is  relatively  small.  The  kink  can  be 
smoothed  by  interpolation  using  the  spectral  data  in  the 
neighborhood  of  k  =  1  rad/m.  A  simple  linear  interpolation 
is  applied  to  B(k)  in  the  range  6 kp  <  k  <  1.2  rad/m  using 
values  of  B(6kp)  and  5(1.2).  No  smoothing  is  applied  if 
6kp  >  1  rad/m.  Examples  of  the  smoothed  spectra  are 
shown  in  Figure  2d  in  log-log  scales  for  wind  seas  and 
mixed  seas  together.  The  smoothing  of  the  spectral  kink 
only  introduced  a  minor  change  in  the  mean  square  slope 


computation  and  does  not  alter  the  main  conclusions 
discussed  in  this  paper. 

3.  Mean  Square  Slope 

[8]  Using  the  wave  number  spectrum  presented  in  the  last 
section,  the  mean  square  slope,  s2,  of  the  ocean  surface  can 
be  calculated.  Because  waves  longer  that  the  wavelength  at 
the  spectral  peak  make  only  negligible  contribution  to  the 
total  mean  square  slope,  the  lower-bound  wave  number  for 
integration  is  set  at  kp.  The  highest  wave  number  of  the 
spectrum  considered  here  is  316  rad/m,  thus  the  mean 
square  slope  computed  in  the  following  covers  almost  the 
fidl  range  of  the  gravity  wave  spectrum. 

[9]  The  results  calculated  for  the  upper-bound  cutoff 
wave  numbers  of  316,  63,  and  21  rad/m  (cutoff  wavelengths 
0.02,  0.1,  and  0.3  m)  for  wind  seas  and  mixed  seas  are 
shown  in  Figure  3.  For  comparison,  field  data  of  the  mean 
square  slopes  obtained  from  analyses  of  sun  glitter  [Cox  and 
Munk,  1954]  and  airborne  Ai„-band  radar  backscattering 
cross  sections  [Jaclison  et  al.,  1992]  are  superimposed  in 
the  figure.  The  sun  glitter  data  are  further  divided  into  two 
groups,  clean  surface  and  artificial  or  man-made  slicks.  The 
slicks  suppress  short  waves,  and  the  cutoff  wave  number 
was  estimated  to  be  about  21  rad/m,  that  is,  waves  shorter 
than  0.3  m  were  suppressed  [ Cox  and  Munk,  1954],  The 
agreement  between  field  data  and  computations  appears  to 
be  reasonable  considering  the  large  data  scatter  due  to 
difficulties  in  the  acquisition  and  analysis  of  such  measure¬ 
ments.  The  sun  glitter  data  in  clean-water  conditions  com¬ 
pare  well  with  the  mean  square  slope  integrated  from  kp  to 
3 1 6  rad/m,  which  does  not  extend  into  the  capillary  regime 
of  the  wave  spectrum.  In  later  theoretical  analyses  of  the 
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Figure  4.  The  (a)  slope  b\  and  (b)  intercept  b0  describing  the  linear  wind  speed  dependence  of  the  mean 
square  slope  integrated  from  kp  to  ks.  The  polynomial  functions  of  b  ]  (ks)  and  b0(ks)  for  wind  seas  and 
mixed  seas  are  given  in  equation  (5)  and  illustrated  as  continuous  curves  in  the  figure.  See  color  version 
of  this  figure  in  the  HTML. 


relationship  between  the  surface  wave  spectrum  and  the 
mean  square  slope  [e.g.,  Plant,  1982;  Phillips,  1985],  it  has 
been  speculated  that  capillary  waves  may  not  contribute  a 
significant  portion  to  tire  ocean  surface  mean  square  slope 
but  a  definitive  proof  remains  unavailable.  On  the  other 
hand,  Wentz  [1976]  evaluated  the  assumptions  used  in  the 
data  processing  of  sun  glitter  by  Cox  and  Munk  [1954],  He 
concluded  that  the  contribution  from  very  steep  surface 
slopes  cannot  be  recovered  from  the  mean  square  slope  data 
derived  from  the  sun  glitter  measurements  as  processed  by 
Cox  and  Munk  [1954].  Thus  the  result  reported  by  Cox  and 
Munk  [1954]  only  serves  to  provide  a  lower  bound  of  the 
ocean  surface  mean  square  slope.  The  analysis  of  interme¬ 
diate-scale  waves  presented  here  seems  to  support  the  view 
point  expressed  by  Wentz  [1976]. 

[io]  Jackson  et  al.  [1992]  derived  the  mean  square  slope 
of  the  ocean  surface  from  measurements  of  the  airborne  Ku- 
band  altimeter  backscattering  data.  Using  their  estimation  of 
the  diffraction  limit  of  0.1  m,  the  mean  square  slope 
integrated  from  kp  to  the  corresponding  wave  number 


mately  linearly  with  wind  speed.  The  relationship  can  be 
written  as 

j2(A>)  =  b\ {ks)U!0  +  b0(ks).  (4) 

The  dependence  on  ks  of  the  mean  square  slope  and  the 
coefficients  of  the  linear  function  of  wind  speed 
dependence  (slope  b\  and  intercept  b0)  is  written  explicitly 
in  equation  (4).  By  varying  C/10  from  2  to  20  m/s  in  steps 
of  1  m/s  and  ks  from  10  to  310  rad/m  in  steps  of  10  rad/m, 
a  dataset  of.?  as  a  function  of  Uw  and  ks  is  created.  Least 
squares  fitting  procedure  is  applied  to  compute  b\{ks )  and 
b0(ks)  for  the  wind  seas  and  mixed  seas  (Figure  4).  The 
background  waves  in  mixed  seas  increase  the  ambient 
level  of  the  ocean  surface  roughness  (reflected  by  bo)  and 
reduce  the  rate  of  increase  of  surface  roughness  with  wind 
speed  (reflected  by  b ]).  The  systematic  variations  of  b\(ks) 
and  b0(ks)  can  be  represented  by  the  following  polynomial 
functions 


6.34  x  lO'11^  -  5.26  x  KrtVc2  +  1.98  x  lO"5^  +  1.72  x  10-3,wind  seas 

3.72  x  IO-"*3  -  3.43  x  10"8/t2  +  1.30  x  lO’5^  +  1.75  x  10“3, mixed  seas 

5.72  x  10_14k3  —  5.53  x  10-"/t4  +  2.07  x  10"SA;3  -  3.79  x  10-6k2  +3.31  x  10 ~Aks  -  1.13  x  10“3,wind  seas 
7.41  x  lO-'4^5  -  7.06  x  10~”k4  +  2.60  x  10”8k3  -  4.67  x  lO"6*2  +  4.22  x  10“%.  -  9.39  x  10^4, mixed  seas 


•  (5) 


(62.8  rad/s)  yields  reasonable  agreement  with  their  analysis  Equation  (5)  can  be  used  to  calculate  the  mean  square 
(Figure  3).  slopes  of  arbitrary  wave  number  bands  of  surface  gravity 

[li]  Empirically,  the  mean  square  slope  integrated  from  waves.  The  quantification  of  the  band-pass  filtered  mean 
kp  to  an  upper  bound  wave  number,  ks,  increases  approxi-  square  slopes,  such  as  the  tilting  and  Bragg  scattering 
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Figure  Al.  (a)  The  dispersion  relation  of  surface  waves  in  the  presence  of  background  currents,  (b)  The 
ratio  between  the  current-induced  Doppler  frequency  shift  and  the  intrinsic  wave  frequency.  Two- 
dimensional  spectra  obtained  by  wave  gauge  arrays  in  (c)  fixed-station  and  (d)  free-drifting 
measurements.  See  color  version  of  this  figure  in  the  HTML. 


components  of  the  sea  surface  roughness,  are  frequently 
needed  in  many  remote  sensing  applications  such  as  the 
computation  of  the  radar  scattering  cross  section  and  the 
modulation  transfer  function  [e.g.,  Plant,  1986,  1997]. 

4.  Summary 

[12]  Field  measurements  of  intermediate-scale  waves 
were  obtained  using  a  free-drifting  measurement  technique 
to  reduce  the  complication  in  frequency-to-wave  number 
conversion  associated  with  Doppler  frequency  shift  encoun¬ 
tered  in  fixed-station  data.  Following  an  analysis  by  Phillips 
[1984]  on  the  source  function  balance  of  surface  waves,  the 
dependence  of  the  dimensionless  spectrum,  B(k),  on  the 
dimensionless  wind  friction  velocity,  u*/c,  is  investigated. 
The  result  shows  that  B(u*/c)  can  be  represented  by  a  power 
law  function.  The  proportionality  coefficient,  A0,  and  the 
exponent,  a0,  of  the  power  law  function  are  determined 
from  field  data.  The  parameterization  functions  A0(k)  and 
a0(k)  can  then  be  used  to  construct  the  wave  number 
spectra  of  intermediate-scale  waves  at  different  wind  speeds 
(section  2).  Using  this  spectral  model,  the  band-pass  filtered 
mean  square  slope  is  derived.  The  calculated  mean  square 
slopes  are  in  good  agreement  with  those  obtained  from 
analyses  of  sun  glitter  [Cox  and  Munk,  1954]  and  Ku 
band  altimeter  scattering  cross  sections  [Jackson  et  al., 
1992],  Empirically,  it  is  found  that  the  mean  square  slope 
integrated  from  kp  to  ks  (k,  >  ~10  rad/m)  increases 
approximately  linearly  with  wind  speed.  For  convenient 
computation  of  the  band-pass  filtered  ocean  surface  rough¬ 


ness  components  frequently  needed  in  remote  sensing 
applications,  polynomial  functions  on  the  wave  number 
dependence  of  the  slope  and  intercept  of  the  linear  wind 
speed  function  are  presented  as  equation  (5). 

Appendix  A:  Doppler  Frequency  Shift 

[13]  Figure  Ala  shows  the  dispersion  relation  w  =  w0  + 
Aw,  and  Figure  Alb  shows  the  ratio  between  the  current- 
induced  Doppler  frequency  shift,  Aw  =  u  •  k,  and  the 
intrinsic  frequency,  w0  =  (gk  +  Tk3)1'2,  where  u  is  the  surface 
current  vector  with  modulus  u,  k  is  tire  wave  number  vector 
with  modulus  k,  and  r  is  the  ratio  between  the  surface 
tension  and  the  water  density.  The  range  of  the  surface 
currents  used  in  the  computation  is  from  -0.5  to  0.5  m/s, 
which  is  a  common  range  of  wave-induced  orbital  velocities 
in  field  or  laboratory  environment.  For  example,  current 
amplitude  of  0.5  m/s  can  be  generated  by  waves  with  6-s 
wave  period  and  0.5-m  wave  amplitude,  or  1-s  period  and 
0.08-m  wave  amplitude.  As  shown  in  Figure  Alb,  the  signal 
from  Doppler  frequency  shift  is  more  than  twice  the 
intrinsic  value  of  gravity-capillary  waves.  The  Doppler 
frequency  shift  is  still  quite  substantial  in  the  meter-long 
wave  components  and  amounts  to  about  40  percent  of  the 
intrinsic  frequency  at  ±0.5  m/s  current  levels.  Hwang  et  al. 
[1996]  suggested  that  the  problem  of  Doppler  frequency 
shift  can  be  considerably  reduced  by  using  a  free-drifting 
measurement  technique.  Experiments  were  carried  out 
to  compare  the  wave  number-frequency  spectra,  S(k,  w), 
measured  by  linear  arrays  of  wire  gauges  deployed  in  free- 
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drifting  and  fixed-station  configurations  in  a  canal  about 
100  m  wide  and  400  m  long.  With  fixed-station  measure¬ 
ments,  the  distribution  of  S(k,  to)  show  the  typical  shift  from 
the  linear  dispersion  relation  (Figure  Ale).  The  magnitude 
of  the  net  frequency  shift  is  in  agreement  with  that  induced 
by  a  surface  drift  current  of  about  two  to  three  percent  wind 
speed.  There  are  also  apparent  upwind-traveling  wave 
components  (negative  k  and  positive  to  in  Figure  Ale).  In 
contrast,  the  distribution  of  S(k,  to)  from  free-drifting  mea¬ 
surements  concentrates  along  the  linear  dispersion  curve  and 
the  upwind  traveling  components  disappear  (Figure  Aid). 

[14]  Acknowledgments.  This  work  is  sponsored  by  the  Office 
of  Naval  Research  (NRL  PE  61153N  and  62435N).  (NRL  contribution 
JA-7330-05-5 1 99.) 

References 

Cox,  C.  S„  and  W.  Munk  (1954),  Statistics  of  the  sea  surface  derived  from 
sun  glitter,  J.  Mar.  Res.,  13,  198-227. 

Hara,  T.,  E.  ,T.  Bock,  J.  B.  Edson,  and  W.  R.  McGillis  (1998),  Observation 
of  short  wind  waves  in  coastal  waters,  J.  Phys.  Oceanogr.,  28,  1425- 
1438. 

Hwang,  P.  A.,  and  D.  W.  Wang  (2004),  An  empirical  investigation  of  source 
term  balance  of  small  scale  surface  waves,  Geophvs.  Res.  Lett.,  31, 
LI 530 1 ,  doi:  1 0. 1 029/2004GL020080. 

Hwang,  P.  A.,  D.  B.  Trizna,  and  J.  Wu  (1993),  Spatial  measurements  of 
short  wind  waves  using  a  scanning  slope  sensor,  Dvn.  Atmos.  Oceans,  20, 
1-23. 


Hwang,  P.  A.,  S.  Atakturk,  M.  Sletten,  and  D.  B.  Trizna  (1996),  A  study  of 
the  wavenumber  spectra  of  short  water  waves  in  the  ocean,  J.  Phys. 
Oceanogr.,  26,  1266-1285. 

Jackson,  R  C.,  W.  T.  Walton,  D.  E.  Hines,  B.  A.  Walter,  and  C.  Y.  Peng 
(1992),  Sea  surface  mean-square  slope  from  RT„-band  backscatter  data, 
J.  Geophvs.  Res.,  97,  11,411-11,427. 

Makin,  V.  K.,  and  V.  N.  Kudryavtsev  (1999),  Coupled  sea  surface- 
atmosphere  model  1.  Wind  over  waves  coupling,  J.  Geophys.  Res., 
104,  7613-7623. 

Makin,  V.  K.,  and  V.  N.  Kudryavtsev  (2002),  Impact  of  dominant  waves  on 
sea  drag,  Boundary  Layer  Meteorol.,  103,  83  -  99. 

Phillips,  O.  M.  (1984),  On  the  response  of  short  ocean  wave  components 
at  a  fixed  wavenumber  to  ocean  cunent  variations,  J.  Plivs.  Oceanogr., 
14,  1425-1433. 

Phillips,  O.  M.  (1985),  Spectral  and  statistical  properties  of  the  equili¬ 
brium  range  in  wind-generated  gravity  waves,  J.  Fluid  Mech.,  156 , 
505-531. 

Plant,  W.  J.  (1982),  A  relationship  between  wind  stress  and  wave  slope, 
J.  Geophys.  Res.,  87,  1961-1967. 

Plant,  W.  J.  (1986),  A  two-scale  model  of  short  wind-generated  waves  and 
scatterometry,  J.  Geophys.  Res.,  91,  10,735-10,749. 

Plant,  W.  J.  (1997),  A  model  for  microwave  Doppler  sea  return  at  high 
incident  angles:  Bragg  scattering  from  bound,  tilted  waves,  J.  Geophys. 
Res.,  102,  21,131—21,146. 

Wentz,  R  J.  (1976),  Cox  and  Munk’s  sea  surface  slope  variance,  J.  Geo¬ 
phys.  Res.,  81,  1607—1608. 


P.  A.  Hwang,  Oceanography  Division,  Naval  Research  Laboratory, 
Building  1009,  Room  B151,  Stennis  Space  Center,  MS  39529,  USA. 
(phwang@nrlssc.navy.mil) 


7  of  7 


